Abstract. The excitation of low frequency plasma waves 
in an expanding anisotropic cosmological model which 
contains a magnetic field frozen into the matter and point- 
ing in the longitudinal direction is discussed. Using the 
exact equations governing finite-amplitude wave propaga- 
tion in hydromagnctic media within the framework of gen- 
eral theory of relativity, we show that a spectrum of mag- 
netized sound waves will be excited and form large scale 
"damped oscillations" on the expanding universe. The 
characteristic frequency of the excited waves is slightly 
^««J . shifted away from the sound frequency and the shift de- 

' pends on the strength of the primordial magnetic field. 

, This magnetic field dependent shift may have an effect on 
i * ' the acoustic peaks of the CMB. 

o : 
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1. Introduction 

At present not much is known about the existence or ab- 
sence of magnetic fields in the early Universe. Several sce- 
narios exist for the generation of magnetic fields (see re- 
view by Grasso & Rubinstein 2001). The implications of 
the primordial magnetic fields for the formation and evo- 
lution of the observed structures have been the subject of 
continuous theoretical investigation. 

The observed magnetic field in the clusters and galax- 
ies is in energy equipartition with the gas and the cosmic 
rays in these systems. The observed magnetic field may be 
primordial in nature, being left over from the early uni- 
verse, or created when the first structures formed. Even 
if we assume the presence of a protogalactic dynamo to 
explain the magnitude of the magnetic field still requires 
a small primordial field. It is a difficult task to determine 
if the origin of the cluster or the galactic magnetic field 
is primordial. In contrast the role of the magnetic field in 
the fluctuations spectrum of the cosmic microwave back- 
ground anisotropics can provide bounds in the large scale 
magnetic field of the early universe (see Adams et al. 1996, 
Barrow et al. 1997, Enqvist 1997, Olinto 1997, Durrer et 
al. 1998). 

Many recent studies used a Newtonian or a Friedmann- 
Robertson- Walker (FRW) model for the evolving Universe 
and super-impose a large scale ordered magnetic field. 
The magnetic field is assumed too weak to destroy the 
Robertson- Walker isotropy. The anisotropy induced by 
the magnetic field is treated as perturbation (Durrer et al. 
(1998), Ruzmakina & Ruzmakin 1971, Tsagas & Barrow 
1998). Current observations give a strong motivation for 
the adoption of a FRW model (see reviews above) but the 
uncertainties on the cosmological standard model are sev- 
eral and the limits of the approximations and the effects 
one losses by neglecting the anisotropy of the background 
magnetic field should be investigated. 

In this article, departing from the traditional ap- 
proach, we adopt an anisotropic cosmological model, de- 
veloped initially by Thorn(1967) and Jacobs(1968). In 
this model the magnetic field is responsible for the ini- 
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tial anisotropic expansion of the universe. The model is 
homogeneous and has two equivalent "transverse" direc- 
tions and one inequivalent "longitudinal" direction at each 
point in space time. It contains a perfect fluid obeying a 
"stiff" equation of state p = 7p, with 1/3 < 7 < 1. The 
magnetic field is frozen into the matter and directed along 
the longitudinal direction. 

We analyze the stability of the linear perturbations 
using the above model. Magnetized density perturbations 
were studied in detail in |l8| using a Newtonian model 
for the evolving Universe. Several authors used the FRW 
model to study magnetized perturbations in the radia- 
tion era (see Holcomb & Tajima 1989), the dissipation of 
MHD waves in recombination era (Jedamzik et al.2000) 
and in the inflanatory era (Brevik & Sandvik 2000). The 
propagation of the magnetosonic waves in sub-horizon and 
super-horizon scales were discussed by Tsagas & Maartens 
(2000a) and more recently a magnetized Biannchi I back- 
ground model was used to study the coupling of the mag- 
netism with the geometry (Tsagas & Maartens 200b). 

The main scope of this study is to search for the spec- 
trum of the unstable low frequency plasma waves in an 
anisotropic model for the universe, using the formalism 
developed by Papadopoulos & Esposito (1982). In section 
2 we present the basic equations used and in section 3 
we perturb the MHD equations keeping first order terms 
and searching for solutions when the perturbed quantities 
have the form of a plane wave. We derived the linear dis- 
persion relation for the magnetosonic waves and search for 
the spectrum of excited waves. 

2. Basic equations 

The general relativistic version of the magnetohydrody- 
namic equations are: 

H 2 

xu a + xii a + x6u a + (p+ —);bg ab 

- {H a H b ), b = 0, (1) 
(e - ^-)- ab u a u b = h ab (p + ^-), ab + 2{H 2 6) 
9ft 2 

- (H a H b ), ab + 2x(^- + a 2 -u; 2 - u a ii a ) 
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(p + 3p + H 2 ) + 2u a {H a H b 



2 

+ (H 2 ), a u a 
and 

pH 2 p 



(2) 
(3) 

u a is the fluid veloc- 
p + pH, pli is the inter- 
is the expansion velocity, 
x =_e+p + H 2 , H the magnetic field and h ab = g ab + u a u b . 



49 pH 2 



8ir 4n 
where c—G—l,p is the permeability. 

ity, p is the mass density, e = 

nal energy, ii a = u a c u c , 9 = u a a 



In H the Eqs. ©,(§), and (|) were perturbed using also 
the the condition 6g a b = known as Cowling approach 
|(| and keeping only first order terms in pressure, density, 
velocity, and magnetic field (see Eqs (47), (48) and (49) 
in Papadopoulos & Esposito 1982). 

We linearize the perturbed and search for the am- 
plification or damping of small amplitude hydromagnetic 
waves in the early universe described by the well known 
anisotropic cosmological model (Thorn 1967), 
ds 2 = -dt 2 + A 2 (dx 2 +dy 2 ) + W 2 dz 2 ; (4) 
here, (i) the model contains perfect fluid obeying the equa- 
tion of state p = 7p where | < 7 < 1, (ii) the fluid is as- 
sumed to be co-moving with the coordinate system and 
therefore it M = (—1,0,0,0) and (iii) as seen in the rest 
frame of the fluid there is a magnetic field of strength 
H pointing in the z-direction but no electric field. The 
functions A and W which enter into line element are 
A = A(t) = t 1 ' 2 and W = W{t) = t\ while 
3-7 



167rf 2 (l + 7)' 
(l- 7 ) 1/2 (3 7 



1 



,1/2 



u 

2i(l + 7 ) 
and I = (l- 7 )/(l + 7). 

The ratio of the cyclotron (ujb) to plasma (u p ) fre- 
quency is a function of 7 (uj 2 /uj 2 ) = (2(1 — 7)(37~ l))/(3 — 
7). The cyclotron frequency tends to zero as 7 — > 1 or 
7 —v 1/3 and the ratio < u 2 /ui 2 < 0.3 remains bounded 
when 1/3 < 7 < 1. 

The Standard Hot Big Bang (SHBB) is roughly di- 
vided into two regimes, the radiatin dominated era, for 
which the 7 = 1/3, and the (pressurless) matter domi- 
nated era, for which 7 ~ 0. Thus, at first, it seems not 
possible for the Universe to experience a stiffer equation 
of state, with 1/3 < 7 < 1. However, is simple extension 
of SHBB it is possible to acchieve this in a rather natural 
way. 

One prominent example is cosmologies which include 
a scalar field component in the matter content of the Uni- 
verse. A homogeneous scalar field <f> = 4>(t) can be treated 
as an ideal gas with density p^ — pkin + V and pressure 
P<t, = Pkin - V, where pkm = (l/2)</> 2 and V = V(4>) are 
the field's kinetic and scalar potential energy respectively. 
In the case V » pkin, then 7^ ~ —1 and, if dominates 
the Universe then the latter is forced to undergo acceler- 
ated expansion. Such a situation, if it taken to occur at 
the very first stages of the Universe evolution, is the typi- 



cal realization of inflation, in which <f> is called the inflaton 
field (Olive 1990). At the end of the inflationary period <fi 
usually oscillates around the minimum of V(4>) and de- 
cays into other particles creating the thermal bath of the 
SHBB. However, there exist non-oscillatory models of in- 
flation in which the inflation does not decay at the end of 
the inflation but, instead, it engages into a rapid roll-down 
of its steep potential (Felder et al. 1999). In such models 
during this phase the Universe continues to be dominated 
by the scalar field but this time we have V << p^in so 
that 70 ~ l.This period is usually refereed to as kination 
or deflation (Joyce & Procopec 1998). Eventually kination 
ends and the SHBB begins. At the transition between ki- 
nation and radiation domination we have a change of 7 
from 1 — > 1/3. 

In the next section we derive the dispersion relation for 
the low frequency waves, assuming that all the perturbed 
quantities (Sp, Sp, SH^, <5m m ) can be expressed in the form 

gi(nt — kz) 



3. Excitation of low frequency plasma waves 

We assume that the wave vector and the external mag- 
netic field are along the z-axis = (0, 0, H 3 ) , fc M = 
(0,0, 0,fc), and the perturbations of the fluid velocity 
and the magnetic field have a general form 8H^ = 
(8H l ,5H 2 , SH 3 ), Su^ = (Su°, Su^.Su 2 , 5u 3 ). After some 
long but straight forward calculations, we derive the gen- 
eralized dispersion relation, 



n 2 + k 2 c 2 s -(l + c 2 s ){-r 3 9 2 + 2<j 2 + ^(p + 3p) 
2H 2 9 A9 2 H 2 2H 2 ^ 



p + p 3(p + p) 
l -{p + p){l + Zc 2 s ) 



in [4(9 



2H 2 



3 



A9H 2 



(p + p) 3{p + P )' 
(Ki+inK 2 )- H 



Se ' 



H 3 



-[1 - Ai - mA 2 ], 



and 
SH 3 

~5~T ~ 4tt(2-7)' 
where the functions TZx, IZ2, A±, A2 and T> are: 



Hi = [#3,oo - 4(ftff 3 ,o + 0H3) - 2ff 3 r^] 
+ H 3 [-^H 2 + ^9 2 + 2a 2 + 2{p + 2p)} 

+ ff 33 [-2r33 - 2r3 3 r^ + 2r 3 l 3 r3 ]if 3 



(5) 
(6) 



(n 2 + k 2 )H 3 , 
iff 
1 



n 2 = 2ff 3 ,o - 4#ff3 - r§ 3 ff35 33 , 



Ai = ^[k 2 (H 3 ) 2 + (2 1 -l)9[2xT 3 + (p + H 2 /2),o] 

k 2 c 2 s9 33 (e+p)}, 
A 2 = (2 7 - l)x0, 
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-n l x + k\H 6 y - (2-7)0 



(P- 



— ),o 



rr2 

— -(2-7)^6 



In the absence of magnetic field (e.g 7 = 1/3) the 
model used here reduces to a the standard FRW model 
without cosmological constant, since A(t) = W(t) = t 1 / 2 
and the real part of the dispersion relation (Eq. (|J)) has 
a very simple form 

(i + «£)[!** 



= k 



2 c 2 



2 yP 



L 3 
P)(l 



•3p)] 



3cf 



(7) 



Eq.(|7j) reduces further if we use an Enstein cosmalogical 
model with cosmological constant A = (l/2)(p+3p), which 
makes the medium static and homogeneous and 9 = 0. Us- 
ing the above assumptions the Jeans instability (n 2 < 0) 
for a static isotropic and homogeneous model was derived 
by in Jackson (1972). The Jeans instability is quenched by 
expansion (9 ^ 0) in the FRW model, without cosmologi- 
cal constant, since the second term in the right hand side 
of Eq. (Q) is always positive and dominates over the third 
which drives the Jeans instability in the static model. 

It is easy to verify that the dispersion relation (Eq. |J) 
is simplified considerably in the limit t — ► 0, 



n 2 \l 



k 



1 



2-7' ' "2-7 
+ T 1 ( 1 ,t)+inT 2 ( 7 ,t) = 0, ; ( (8) 

where c s is the sound speed and u 2 A = H 2 /Airp is the 
Alfven speed. The functions T\ and T 2 have the following- 
functional form 

T l(7 ,t)- fl(7) 



f 2 



= -(9 7 6 - 993 7 5 + 4842 7 4 - 68I87 3 + 2889 7 2 
- 461 7 + 788)/(6( 7 - 2)( 7 - 3) 2 (1 + 7 ) 2 ) 



and 

r 2 (7,*) 



-40+ 4 1( f , 
3 \4ir(p + p) 



A9 



H 2 



4tt(2 - 7) 



2H 2 



Air(p+p) 
(1-7) ( 



4?rp(2 
H 2 



7) 



(9) 

(7+l)t \^p{2- 1 )) ' v ; 

where 9 = — 2/((l + 7)£).The function T2 depends on 7 
and t since the magnetic field, the pressure and density 
are all functions of 7 and t. In other words Ta(7, t) has 
the form 

where 
C(7) = 



t 



457 4 - 367 3 - I267 2 + 212 7 - 127 



3( 7 -3)( 7 -2)(7+l)2 
We can easily show that Ti(7, t) is positive and Tzfo, t) 
negative for all allowed values of 1 /3 < 7 < 1. 

We re-write Eq. (||) in the form D r + iDi = and 
assuming that the the real part of the excited frequency 
is much larger than the imaginary part, n r >> n il we 



obtain the frequency of the excited wave from the equation 

Dr = 0, 



nl\l 



7 



-] = k 



2 c 2 



2-7 



t 2 



(10) 

and the imaginary frequency from the well known rela- 
tion m = —Di(n r , k) /(dD r (n r , k)/dn r )(see Krall & Triv- 
ielpiece 1973) 

Ti ■ : 

4 2(1-^(^/(2-7)) 

[Ch)\ ( 1 

\ t J Ul-t&(*)/(2-7)) 
were has a large negative value when t 



(11) 



and decay 

fast at later times. The role of magnetic field is to make the 
expanding universe even more stable against the expected 
Jeans type instabilities discussed early. 

The real frequency on the other hand is shifted away 
from the sound speed 

+ 



k 2 c 2 ' ' 



2-7 



11 



1 



(12) 



as 7 — -> 1/3 the magnetic field approaches zero and the 
anisotropic model used here approaches the weakly mag- 
netized FRW model used extensively in the literature. The 
characteristic frequency n r will approach the frequency 
n r = kc s + An r (uA)- The presence of the Alfven speed in 
the frequency n r in Eq. (|i~2| ) may be responsible for the 
distortion of acoustic peaks, as it was pointed out first by 
Adams et al. (1996). 

The imaginary part is negative rij < and a spectrum 
of low frequencies will be excited with frequency n r . The 
linear density perturbation will have the form 

^~[p e |5h '' )l |e i ^ t - fa ' (13) 

wereB( 7 ,i) = C( 7 )/(2(l - u 2 A (t)/(2 - 7))). |B( 7 ,t)| is a 
decaying function of time. For t — ► 0, 7 — > 1 then ua —* 1 
and \B(j,t)\ — > 00. As t >> and 7 — > 1/3 then ua —> 
and |B(l/3,t)| — > (3/2). We then conclude that in the 
early universe the amplitude becomes large and the ex- 
cited waves will form a spectrum of damped oscil- 
lations in the expanding magnetized universe. We can 
search for the source of damping by analyzing the func- 
tional form of T2( 7 , t) in Eq. (||). In the absence of mag- 
netic field (7 = 1/3) T2 = —49 and the expanding plasma 
stabilizes the ion acoustic waves, on the other hand if the 
universe is static (A = constant and 9 ~ A/ A ~ 0) the 
stabilization of the magnetosonic waves is caused by the 
magnetic pressure. In a magnetized expanding universe 
(Eq. (0)) both processes will be combined to damp the ex- 
cited large amplitude oscillations. Assuming, for the sake 
of argument that the universe at a certain time i=const is 
becoming static, a Jeans type instability sets in so we can 
easily realize that the damping is caused by the free en- 
ergy available both in the plasma and magnetic pressure 
during the initial stages of its expansion. 
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4. Summary 

We analyzed the wave propagation in the early Universe 
using an anisotropic cosmological model which contains 
a magnetic field frozen into the perfect fluid with equa- 
tion of state p = 7/5 (with 1/3 < 7 < 1) and pointing 
in the z-direction. The choice of a stiff equation of state 
with 1/3 < 7 < 1 may be realistic representation for the 
evolving universe in the phase between the kination and 
the radiation era. 

The anisotropic model for the early universe used 
can approach the weakly magnetized FRW model when 
7 — > 1/3. In this case the anisotropic model used here ap- 
proaches the isotropic FRW model were the magnetic field 
is a small linear perturbation (Durrer et al.1998, Ruzmak- 
ina & Ruzmakin 1971). Our main results in this article are: 
(1) The Jeans instability, found early for an Einstein uni- 
verse with very particular cosmological constant (Jackson 
1972), will be absent since the expansion and the presence 
of magnetic field will act as a stabilizing force. (2) The 
anisotropic magnetized Universe re-enforces the stability 
found initially in the FRW model. (3) A spectrum of low 
frequency and large amplitude damped oscillations will 
appear in the early universe with characteristic frequency 
n r = kc s + An, where An depends on the strength of the 
primordial magnetic field will be excited. This magnetic 
field dependent shift may cause measurable distortion of 
accoustic peaks (see also Adams et al.1996, Jedamzik et 
al. 2000, Koh & Lee 2000, Durrer et al. 2001). The am- 
plitude of the excited waves is large in the early universe 
and gradually decay. 

We then conclude that if the universe have passed 
through a strongly magnetized anisotropic phase, before 
the recombination era started, the waves predicted here 
will be responsible for the formation of large scale fluctu- 
ations, as it has been shown in the 2-D numerical simula- 
tions (Brandenburg et al. 1996). 
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